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ABSTRACT 


The purpose of this paper is to present a method for 
determining the simultaneous solution of an overdetermined 
system of nonlinear equations in R n . A numerical algorithm 
is developed, and some examples of convergent solutions are 
contained in the last section of the paper. The algorithm 
can be used for problems such as obtaining the intersection 
mach cones which is used in the method of characteristics. 

In the development of the solution, it is assumed that each 
function can be expanded in an n dimensional Taylor’s 
Series in a neighborhood of the solution. The method utilizes 
an iterative procedure, involving the inversion of an 
n x n matrix, which does not necessarily converge. Thus, 
the method might be time consuming or inapplicable to 
some problems. 


INTRODUCTION 

The purpose of this paper is to present a method for 
determining the simultaneous solution of an overdetermined 
system of nonlinear functions in R n . Consider a set of 
functions f^:R n ** R n such that, for some y e R n , 
f^y) s 0 for i = 1,2,...,K . Expressing the coordinates in 
R n as (xj ,x 2 , . . . x n ) , it will be assumed that each of the 
functions f i can be expanded in a convergent n-dimensional 
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Taylor's series about any point In some neighborhood of 
y . Being similar to Newton's method for the solution of a 
single equation In one variable, the convergence of the 
solution will be subject to the initial guess for a solution 
and the characteristics of the functions. 

THE SPECIAL CASE OF n =* k 

For the case n ■ k consider the usual approach to this 
problem. Choose a point yo fs a guess to the solution. 
Then, provided all of the Taylor's series converge, find 
Ay such that 

0 * 


0 = 

3f 

+ • • • 3x^ Ax n + • • • 


r x (y o + Ay) - fi(yo) + 


3 f i ( y o ) 

3xi 


Ax 


3fl 

3x? 


Ax? + . . . + 


3fl 

3x“ Ax n + 
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3 f 3 f 

f n (y o + A *> + TT Ax > + 3x7 AX2 
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If the notation 

3 Myo) 

A(yo) * {a ij } 3 ax 1 * 1 » n » J * 1 » n » and 

J 

Ax^ 1 ^ * (Ax*, Ax 2 ,...Ax n ) is used then a first order 

approximation for Ax^ 1 ^ Is given by 


Ax 


(l) 


-l 


(yo)i 


■f i (yo) 
Cy o > 


-f n (yo) 






In the usual fashion an iterative procedure can then be 

( o_ i ) 

initiated by setting y^ * y^ + Ax ' . 

Because each of the series was truncated, the value of Ax 
is an approximation and the method might not converge. 

THE OVERDETERMINED CASE 

For the case k > n , a solution for the equation 
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does not necessarily exist, since A is a k * n matrix. 

Since Ax^ 1 ^ is an approximation even for the case n * k , 

consider the solution Ax^ ' such that £ (p. - (-f.(yo)) 2 

i-i 1 1 

is minimized. If a unique solution exists, then it is given 

by 


'-fiCyoV 


Ax 


( 1 ) 


(A 


(yo) 


A (yo) ) 


-1 


(y 



(yo), 


An iterative solution can now be introduced Just as before 
with 

h * + AxU_1) • 


Mr. Melvin J. Arldt of LEC has written a computer 
program using this method. Development and testing of 
this program has led to the following examples. 

All execution times are given for the UNIVAC 1108. 
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Example 1. 

2 2 

fi ■ X 2 + Xi - 1 

f 2 * X 2 - X! ♦ 1 

f 3 » X 2 ♦ Xi - 1 


Initial guess 

Xi * 

1/2 


x 2 » 

' 2 

Converged in 7 iterations 


Computer time required 

.231 

sec . 

Example 2. 



1*1 * 

2 ' 
x 2 + 

2 

x i — 9 

f 2 * 

xi - 

3 

f*3 * 

X 2 - 

xi + 3 

f4 » 

2 

X 2 + 

(xx - 6) 2 - 9 

Initial guess 

x; 

= 1 


X 2 

* 0 


Converged In 6 Iterations 
Computer time required .22*4 sec. 
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Example 3. 


2 2 2 


fl 

« 

X 1 

«• 

x 2 ♦ X3 - 1 

fj 

■ 

2 

Xj 

■f 

x 2 ♦ (x 3 - 2) 2 - 1 

fj 

M 

Xl 

+ 

x 2 + x 3 - 1 

f. 

■ 

Xl 

♦ 

X 2 « X 3 ♦ 1 

f* 

■ 

3 

Xl 


3 x 2 ♦ ( 5 x 3 - xi ♦ l) 2 


Initial guess xi - 1 

X2 * 2 

X 3 « 1 

Converged in 25 iterations 
Computer time required 1.131 sec. 

The remaining examples use the recursion formula 

fl 

i ■ 2, . . .n - 1 ; f, 

f n 

with a = -.5 8*1 


-(3 + ax! ) x i + 2 x 2 - 8 


x i-i " (3 + ax t ) + 2x 1+1 


x n-i - (3 + ax n } x n ’ 6 
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Example . 


n ■ 5 


Initial guess x 1 ■ -1 i ■ 1,2,.. .5 

Converged in 5 iterations 
Computer time required .508 sec. 

Example 5. 


n - 10 

Initial guess x i * -1 i ■ 1,2,.. .,10 

Converged in 5 iterations 
Computer time required 1.799 sec. 

Example 6. 


n * 20 

Initial guess x 1 * -1 i * 1,2,. ..,20 

Converged in 6 iterations 
Computer time required 10.709 sec. 
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